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Abstract We discuss the temporal variations and frequency distributions of solar wind and 
interplanetary magnetic field parameters during the solar minimum of 2007 - 2009 from 
measurements returned by the IMPACT and PLASTIC instmments on STEREO-A. We find 
that the density and total field strength were significantly weaker than in the previous mini- 
mum. The Alfven Mach number was higher than typical. This reflects the weakness of mag- 
netohydrodynamic (MHD) forces, and has a direct effect on the solar wind-magnetosphere 
interactions. We then discuss two major aspects that this weak solar activity had on the mag- 
netosphere, using data from Wind and ground-based observations: i) the dayside contribu- 
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tion to the cross-polar cap potential (CPCP), and ii) the shapes of the magnetopause and bow 
shock. For i) we find a low interplanetary electric field of 1 .3 ± 0.9 mV m -1 and a CPCP of 
37.3 ± 20.2 kV. The auroral activity is closely correlated to the prevalent stream-stream in- 
teractions. We suggest that the Alfven wave trains in the fast streams and Kelvin-Helmholtz 
instability were the predominant agents mediating the transfer of solar wind momentum and 
energy to the magnetosphere during this three-year period. For ii) we determine 328 magne- 
topause and 271 bow shock crossings made by Geotail , Cluster 1, and the THEMIS B and 
C spacecraft during a three-month interval when the daily averages of the magnetic and ki- 
netic energy densities attained their lowest value during the three years under survey. We use 
the same numerical approach as in Fairfield’s (/. Geophys. Res. 76 , 7600, 1971) empirical 
model and compare our findings with three magnetopause models. The stand-off distance of 
the subsolar magnetopause and bow shock were 11.8 Re and 14.35 Re, respectively. When 
comparing with Fairfield’s (1971) classic result, we find that the subsolar magnetosheath 
is thinner by ~1 R E . This is mainly due to the low dynamic pressure which results in a 
sunward shift of the magnetopause. The magnetopause is more flared than in Fairfield’s 
model. By contrast the bow shock is less flared, and the latter is the result of weaker MHD 
forces. 


1. Introduction 

The exceptionally long and pronounced minimum in solar activity at the transition from So- 
lar Cycles 23 to 24, lasting from 2007 to 2009 inclusive, has excited a lot of attention in the 
community (see, for example, the compilation of articles in Cranmer, Hoeksema, and Kohl, 
2010). Comparisons and contrasts have been drawn with other historic cases of protracted 
low solar activity, as judged for example by the sunspot number, such as the Dalton mini- 
mum (1790- 1830). It has also been noted that unusual solar behavior may have started well 
before the 2007 - 2009 minimum (Feynman and Ruzmaikin, 201 1). It was hypothesized that 
the years 2007 - 2009 represented a minimum of the Centennial Gleissberg Cycle, the last 
minimum of which occurred at the beginning of the 20th century, rather than a Maunder-type 
minimum (Feynman and Ruzmaikin, 2011). 

Diminished solar activity will affect the properties of the solar wind and hence the types 
of interaction it has with planetary magnetospheres. In this paper we concentrate on the 
variation of the solar wind parameters and, in particular, of the energy budget in the inner 
heliosphere (kinetic and magnetic energy densities). Examining their effects on the terres- 
trial magnetosphere, we focus on two key aspects: 

i) an estimate of the cross-polar cap potential in its dependence on these interplanetary 

parameters (, i.e . the dayside source of plasma convection in the magnetosphere), and 

ii) the shapes of the bow shock and magnetopause. 

We also compare these findings to others reported for more typical solar activity, such 
as the minimum in 1995. Results on aspect i) will allow an educated guess as to the solar 
wind-magnetosphere interaction processes most common during this three-year period. Re- 
sults on ii) - admittedly restricted to a three-month interval - will allow us to see in what 
way the magnetopause boundary and its shape (i.e. stand-off distance and flaring) reflect 
the interplanetary parameters and how they differ from three widely used models of the 
magnetopause. 

For considerations of energies in the interplanetary medium we discuss field and plasma 
measurements made by STEREO-A (ST- A). For the impact on the terrestrial magnetosphere 
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we use magnetic field and plasma data from the Wind spacecraft. In 2007 - 2009 Wind was 
orbiting around the Lagrangian LI point ~220 Re upstream of Earth. ST-A, launched in 
October 2006, is in an Earth-like orbit but leads the Earth, increasing its longitudinal sepa- 
ration from Earth by 22.5 “/year. At the end of 2009 its angular separation from Earth was 
~ 68°. The magnetospheric response is monitored by data from Cluster 1, Geotail, and the 
THEMIS B and C spacecraft. 


2. Solar Wind Observations: STEREO-A 

We look at solar wind observations from the vantage point of ST-A. Figures 1-5 give an 
overview of key solar wind and IMF parameters. The magnetic field and plasma data are, 
respectively, from the IMPACT (Luhmann et al., 2008) and PLASTIC (Galvin et al, 2008) 
instrument suites. The format of the figures is the same. The top panel shows the tempo- 
ral profile of the relevant quantity ( B , V, T, N , B N ) at 1 -min resolution. The data set is 
subdivided into three samples (S1-S3) as indicated. SI: March 2007 - December 2007; S2: 
January 2008 -October 2008; S3: November 2008 -October 2009. The subdivision into 
three intervals of approximately equal duration is arbitrary and is merely meant to bring 
out any evolutionary trends that may have occurred. The second panel shows the frequency 
distribution of the relevant quantity over each sample in the form of histograms. Finally, 
the third panel shows a spectrogram of the distribution, arranged by Carrington Rotation 
(CR) and month/year. The occurrence frequency is color-coded as shown on the right of this 
panel. 

Figure 1 refers to the magnetic held strength. All three distributions are skewed to the left 
and exhibit long tails. The distribution is thus not Gaussian and can reasonably well be fitted 
by a log-normal distribution (not shown). The held peaks at 3.5 nT and the spectrogram 
shows that the most probable values lie in the range ~ 2.5 - ~4.5 nT. The distribution is 
narrowest in S2 and the peak is slightly shifted to higher values there. We determine the 
“most probable values” for the physical parameters plotted in Figures 1-8 from a visual 
inspection of the bottom panels of the respective hgures. 

We next consider plasma (proton) parameters. Figure 2, 3, 4 show results for the bulk 
speed ( V), temperature ( T ), and density (A), respectively. The time prohles and histograms 
for V indicate multiple peaks in SI and S2, also rehected in the wide spread of values in 
the spectrogram. In fact, the time prohles for SI and S2, too, are indicative of a succession 
of alternating slow and fast streams of large amplitude with associated stream-stream inter- 
actions (co-rotating or stream interaction regions, CIRs/SIRs; Gosling and Pizzo, 1999), a 
pattern which stands out most clearly in S2. By contrast, the velocity prohle of S3 is more 
spiky and the velocity range is more restricted, with few cases of bulk speeds exceeding 
600 kms -1 . There is a much thinner, two-peak structure in S3, centered on lower values 
(slow solar wind). It may be noted that these two V-peaks in S3 are very similar to the two 
peaks observed at low velocities in the distributions SI and S2 - one below 300 km s -1 and 
the other above 300 kms -1 . So, two peaks at low velocities seem to persist throughout the 
whole period. We thus expect to see also several stream-stream interactions where a slow 
stream (dehned here as V < 450 kms -1 ) overtakes an even slower stream. The multiple 
peaks at low V likely suggest different sources of the slow solar wind. In addition, SI and 
S2 have a high-velocity component with extra peaks. 

The proton temperature distributions (Figure 3) are skewed and possess long tails. Sam- 
ples SI and S2 indicate a broad maximum, possibly a two-peak feature. The distribution in 
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STEREO— A : Magnetic field strength 
1 —minute Impact data 
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Figure 1 The distribution of the total interplanetary magnetic field at STEREO-A. The top panel shows the 
time series, the middle panel gives its distribution over the three times intervals (samples); and the bottom 
panel shows a spectrogram of its distribution as a function of Carrington rotation and month/year. 


S3 is shifted to lower values. Most probable values lie in the ranges 2.5-12 x 10 4 K (SI), 
3-16 x 10 4 K (S2), and 2-9 x 10 4 K (S3). 

Figure 4 shows results for the proton density. The distributions are strongly skewed to the 
left and maxima are reached at very low values, particularly on SI and S2. The narrowest 
distribution occurs in S2. The highest densities are reached in S3 because, as we see from 
the distributions of each sample, the most probable value is shifted to higher values from SI 
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STEREO— A : Velocity 
1— minute Plastic data 
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Figure 2 Similar to Figure 1, but for the solar wind bulk speed. 


to S3 (bottom panel), although the differences between SI and S3 are relatively small. Most 
probable values lie in the ranges 1-5 cm -3 (SI); 1-3 cm -3 (S2), and 2-6 cm -3 (S3). 

Figure 5 shows the distribution of the B N component of the magnetic held. RTN coordi- 
nates are used. In these coordinates unit vector R points from the Sun to the spacecraft (here, 
ST- A), vector T points in the direction of planetary motion, and unit vector N completes the 
right-handed RTN system and is such that the RN plane contains the rotation axis of the 
Sun. The distributions in all three samples are similar and may be modeled by Gaussian 
(normal) distributions centered around 0. This is quite typical behavior. Not so typical are 
the narrow profiles (compare with Figure 13 in Luhmann et al. (1993) for solar minimum 
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STEREO— A : Temperature 
1— minute Plastic data 
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Velocity distributions over Carrington rotation 2053-2089 




Apr Jun Aug Oct Dec Feb Apr Jun Aug Oct Dec Feb Apr Jun Aug Oct 
2007 2007 2007 2007 2007 2008 2008 2008 2008 2008 2008 2009 2009 2009 2009 2009 


Figure 3 Similar to Figure 1, but for the solar wind proton temperature. 


conditions in 1975), as may also be seen from the thin slice around 0 nT where the most fre- 
quent values are clustered and which is of width ~ 3 nT (~ — 1.5 to 1.5 nT). This behavior 
is similar for GSM B : at Wind (not shown). Since B N ( B z ) is the major component in deter- 
mining the magnetic shear at the magnetopause, and hence the likelihood of reconnection 
there, the ST-A and Wind data suggest that reconnection processes, while certainly present, 
were not a dominant aspect of the interaction of this solar wind with the magnetosphere. 

We support the last statement further by a study of solar wind-magnetosphere coupling 
functions. To be specific, we consider the Kan-Lee reconnection electric field (E KL , Kan and 
Lee, 1979; see also Sonnerup, 1974; Hairston, Hill, and Heelis, 2003), which includes B z , 
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STEREO— A : density 
1 —minute Plastic data 
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Figure 4 The distribution of the solar wind proton density from ST-A. 


B y and the IMF clock angle, 0 (, i.e ., the polar angle in the GSM Y Z plane). It is defined as 
Fkl = V B t sin 2 (0/2), where B r = (B 2 + B:) l ^ 2 .Weuse OMNI data (King and Papitashvili, 
2005). Quantity E KL follows a log-normal distribution. Over the three-year period 2007 - 
2009 we find that the most probable value Fkl = 0.28 mV/m with a multiplicative standard 
deviation of 0.64. For comparison we did the same analysis for three years of solar activity 
maximum 2001-2003 and obtained 0.61 mV/m with a multiplicative standard deviation 
of 0.66. While the width of the E K L distributions is thus quite similar (multiplicative factor 
0.66 and 0.64), the most probable value in 2007-2009 is about one-half its value in 2001 - 
2003. 
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STEREO-A : Magnetic field — B n component (RTN) 
1— minute Impact data 
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Figure 5 The distribution of the By component of the magnetic held. RTN coordinates are used. 


2.1. Comparison with Solar Minimum 1995 - 1996 

We now compare some aspects of the distributions in 2007 - 2009 with those in the previous 
minimum (Cycle 22/23), from mid 1995 - mid 1996. We start with the three plasma parame- 
ters for which we compare with results given by Burlaga and Lazarus (2000). In 1995 many 
co-rotating interaction regions were present (see for example their Figure 2), very similar 
to the velocity profile in 2007-2008. All proton plasma parameters V , N and T in 1995 
had double -peaked distributions. On this we have agreement in the V and, though less pro- 
nounced, in the T distributions, but this feature is not so evident in the A-distributions for 
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STEREO-A : Magnetic energy density 
1— minute Impact data 
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Figure 6 The distribution of solar wind magnetic energy densities. 


reasons not clear to us. We shall thus compare their most probable values in the two peaks 
for V and in only the major peak for T and N. For 1995, Burlaga and Lazarus (2000) find 
most probable values of 

i) 377 kms -1 and 618 kms -1 for V; 

ii) 5.5 x 10 4 K for T; and 

iii) 8.3 cm -3 for N. 

Comparing with our results, we have similar values for V and T, particularly in SI 
and S2. The large divergence occurs in N where we record much lower densities (~ 3-4 
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STEREO— A : kinetic energy density 
1— minute Plastic data 
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Figure 7 Similar to Figure 6 but for the solar wind kinetic energy densities. 


vs. 8 cm -3 ). The density is thus the interplanetary plasma parameter which deviates most 
from the previous minimum. 

We now consider the magnetic field in the previous solar minimum. For this we use 
1 -min averages from the OMNI database for May 1995 -May 1996. We obtain (average 
and standard deviation) B = 5.4 ± 2.5 nT and for B z : —0.4 ±3.0 nT. For comparison, the 
most probable values of B during the 2007-2009 solar minimum lie in the range [2.5- 
4.5] nT and the those of B z in the range [—1.5, 1.5 nT]. So the magnetic field strength is 
significantly lower in 2007 - 2009. Low magnetic field strengths favor the growth of the 
Kelvin-Helmholtz instability at the magnetopause (see Section 5). 
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STEREO— A : quasi invariant 

1— minute combined Plastic and Impact data 
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Figure 8 Distribution of the ratio of the magnetic-to-kinetic energy densities, i.e., the solar wind quasi-in- 
variant, QI. Note the significant rise during CRs 2285-2287. 


2.2. Solar Wind Energy Densities 

We next discuss the solar wind energy densities and their ratio, the so-called “solar wind 
quasi-invariant”, which was introduced by Osherovich, Fainberg, and Stone (1999). Fig- 
ures 6 and 7 show the magnetic ( E mag = B 2 /2fio) and kinetic ( F kl „ = l/2m p NV 2 ) energy 
densities, respectively, where m p is the proton mass. The distributions of E ma g feature long 
tails and are of the log-normal variety (Feitner et cil., 2009). The most probable values lie in 
the ranges ~ 2 - 6 x 10 -3 nPa(Sl);3-7 x 10 -3 nPa (S2), and 2-7 x 10 -3 nPa (S3). 
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The kinetic energy density distributions in Figure 7 also have long tails and are skewed 
to the left similar to other plasma parameters. The first two data samples feature broader 
distributions, reflecting a similar feature of the velocity distribution (Figure 2). The most 
probable values lie in the ranges 0.3 -0.7 nPa (SI), 0.2-0. 8 nPa (S2), 0.2 -0.6 nPa (S3). 

The ratio E ]m „J E kin - the so-called quasi-invariant, QI - is plotted in Figure 8. The rel- 
evance of quantity QI to our considerations is that it has been shown to correlate very well 
with solar activity as measured by the sunspot number (SSN), better than the individual pa- 
rameters making up this ratio. It is thus a good proxy for solar activity which, furthermore, 
has the distinct advantage of being based on local in situ interplanetary measurements. The 
correlation with SSN has been found to be very good at 1 AU (Osherovich, Fainberg, and 
Stone, 1999), at large AUs in the outer heliosphere (Fainberg and Osherovich, 2002) and, us- 
ing Pioneer Venus (Fainberg, Osherovich, and Stone, 2001) and Helios 1 and 2 data (Leitner 
et al, 2005) also within 1 AU. QI is related to the Alfven Mach number, M A , by QI = M^ 2 . 
QI distributions for the normal solar wind have been modeled well by log-normal distribu- 
tions (Leitner, Farrugia, and Voros (201 1) and references therein). Most frequent values lie 
in the range [0.004, 0.010], corresponding to the M A range [15.8, 10.0], which is typically 
higher than M A ’s in the normal solar wind at 1 AU (see e.g. Slavin and Holzer, 1981). De- 
partures from this trend may be seen in the decrease during CRs 2075 - 2084 and the sharp 
rise in 2085 - 2087, which could reflect some ongoing processes on the Sun (see further in 
Section 3). 


3. Solar Observations 

Figure 9 shows the daily international relative sunspot number (2007 - 2009) and its fre- 
quency distribution. The data are taken from the NASA/Goddard Space Flight Center web- 
site and are computed as a weighted average of measurements made from a network of 
more than 25 observatories (Clette et al., 2007). The deep minimum of Solar Cycle 23 is 
reflected in the low sunspot number. Over the whole period there were very few sunspots, 
or no sunspots at all. 

Figure 10 shows the evolution of the full-disk net magnetic flux density (in pT) of the 
Sun, taken from the Stanford Solar Observatory (see Scherrer et al, 1977), together with its 
frequency distribution. The full-disk integrated data reflect the combined contributions from 
the quiet Sun, active regions and coronal holes, and are a measure of the total “open” mag- 
netic flux on the Sun. The open magnetic flux originates in active regions but in the course of 
a solar cycle it is redistributed over the solar surface by transport processes, eventually form- 
ing the polar coronal holes (Wang, 2009). Coronal holes are low-density regions of the solar 
corona which contain open magnetic fields along which plasma escapes into the interplan- 
etary space forming the high-speed solar wind. The total open magnetic flux varies roughly 
as the Sun’s total dipole strength, which tends to peak a few years after sunspot maximum 
(Wang, 2009). The quasi-periodic short-term fluctuations of the net magnetic flux density in 
Figure 10 are due to the rigid solar rotation of coronal holes with a mean synodic period of 
about 27 days. Due to the slow evolution and long lifetimes (over several solar rotations) of 
coronal holes, the open flux contribution of individual large coronal holes is well reflected 
in the full-disk integrated data. 

Wang (2009), measuring the magnetic flux at 1 AU for which it is reduced by a factor 
of 10 4 , notes that the open flux density in the period 2008 - 2009 (with three-month averages 
of ~ 1.5 nT at 1 AU) is significantly smaller than in previous sunspot minima (~ 2 nT). The 
much smaller fluctuations in S3 (2009) with an almost Gaussian distribution as compared 
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sunspot number — daily averages 
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Figure 9 Evolution and distribution of the daily international sunspot number for the years 2007 - 2009. 


to S 1 and S2 (2007 - 2008) suggest that this trend becomes even more pronounced in the 
year 2009. 

It may be noted that Figure 8 shows a distinct increase in the energy ratio during CR 
2085-2087. This might reflect processes on the Sun. In this regard we note that during 
CR 2076-2079 a coronal hole (CH) which could be attributed to the new Solar Cycle 24 
emerged and subsequently merged with the near northern polar CH. The new-cycle CH 
at high latitudes contributes to a quite slow solar wind speed (maximum speeds of about 
400 kms -1 ). During CR 2085-2087 two new-cycle CHs emerged and merged with the 
polar CH to the north. In addition, a CH with opposite polarity appeared close to this entity. 
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Sun : magnetic flux (daily averages) 
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Figure 10 Similar to Figure 9, but for daily averages of the full-disk net magnetic flux density. 

The complex system of northern polar CH extension resulted in solar wind streams with 
maximum speeds of about 500 kms -1 (Wang et al., 2010). 

4. Major Effects on the Magnetosphere 

4.1 . The Average Cross-Polar Cap Potential in 2007 - 2009 

We now discuss two major effects this solar activity minimum had on the magnetosphere. 
We consider first the solar wind electric field, defined by £sw = |Vsw x Bswl and shown in 
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OMNI Solar Wind Electric Field 2007-2009 
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Figure 11 Over the three-year period, the top panel shows the dawn-dusk motional electric field in the solar 
wind (V X {B r + /> y ) 1 " ) , the second panel shows the cross-polar cap potential (CPCP) after Boyle, Reiff, 
and Hairston (1997); the third panel shows the bulk speed profile in 2008; and the fourth panel auroral, the 
electrojet AL index. Red guidelines connect stream interaction regions with enhancements in AL activity. 


the top panel of Figure 1 1 . We assume that the velocity of the solar wind is in the — x direc- 
tion. One-min averages are used. Average and standard deviation are E = 1.3 ± 0.9 mV/m. 
This is much below typical values in interplanetary coronal mass ejections (ICMEs), which 
are in the range of several mV m -1 to a few tens mV m -1 (e.g. Flairston, Hill, and Heelis, 
2003; Hairston, Drake, and Skoug, 2005). The values of £sw shown in the figure are also 
generally below values when the magnetospheric response to interplanetary driving becomes 
non-linear, which are estimated to be in the range 6-10 mVm" 1 (Siscoe et al., 2002; 
Hairston, Hill, and Heelis, 2003; Muhlbachler et al., 2005). Thus saturation levels in 
magnetospheric response (for example, of the cross-polar cap potential, CPCP) are not 
reached. 

This linear response justifies the use of an empirical formula for calculating the CPCP 
from interplanetary parameters. The CPCP is, of course, a major quantity governing the 
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dynamics of the magnetosphere since it drives the plasma convection cycle. This empirical 
relation was obtained by Boyle, Reiff, and Hairston (1997) from a statistical analysis of 
several years of observations made by the DMSP satellites. It is given by 

<p A = 10“ 4 V 2 + 11.75 sin 3 (0/2) (kV), 

where V is in kms -1 , B is in nT and angle 9 = arccos(B-/R) with B z in GSM coordi- 
nates (i.e., the clock angle of the IMF). This potential gives a contribution to the CPCP, 
which is dependent on the solar wind parameters, i.e. the day side source (as opposed to 
the nightside contribution, which is an independent source of convection and can give a 
significant contribution to the total CPCP during substorm expansion and recovery phases, 
with CPCP increases by a few tens of kV (Lockwood et al., 2009; Kullen et al., 2010; 
Gordeev et al., 2011). The first term depends only on the velocity of the solar wind, while 
the second term is based on the interplanetary magnetic field and its orientation in the GSM 
Y Z plane. 

The time sequence of the “Boyle” potential O a is plotted in Figure 11, panel 2. The 
average value is <t> B = 37.3 ± 20.2 kV. At this high temporal resolution, isolated values 
approach ~ 240 kV, which is typical of saturated CPCP (dependent on solar wind dynamic 
pressure and ionospheric conductivity; e.g. Siscoe et al., 2002; Hairston, Hill, and Heelis, 
2003; Hairston, Drake, and Skoug, 2005). For comparison purposes, typical CPCP values 
lie in the range [40, ~ 80] kV, and a CPCP < 40 kV would correspond to quiet conditions 
[Marc Hairston, private communication, 2012]. Thus the bulk of our values lie in the lower 
half of this range. 

The third panel of Figure 1 1 shows the bulk flows in 2008. As noted earlier, a succession 
of slow and fast streams is evident. In the bottom panel we plot the auroral AL electrojet 
index. One can see that there is practically a one-to-one correspondence between AL in- 
tensifications and CIRs/SIRs (see vertical red guidelines). This suggests that the part of the 
magnetosphere driving due to reconnection comes from the Alfvenic fluctuations in the fast 
streams of stream-stream interaction regions where B z fluctuates about 0 nT (Tsurutani and 
Gonzalez, 1987). 

We now give two examples in our data set of these Alfvenic fluctuations. Figure 12 
shows the first example. For 29 September -8 October 2008 the panels give the proton 
density, temperature, the magnetic field strength and components, the velocity components 
and the auroral AE index. GSE coordinates are used. A very fast stream, reaching above 
700 kms -1 , is overtaking a slower stream. The bottom panel shows the auroral activity. 
We analyze further the interval 1-8 October, marked in the bottom panel by a horizontal 
bar. We form 13 -point averages to represent the background magnetic field and velocity 
vectors. Subtracting this from the measurements gives us the perturbations. From these we 
determine the fluctuations perpendicular to the background magnetic field (ABj_, AVjJ and 
check the Alfven relation component-wise AB ± = (/xop)^ 2 AVx = AA ± ), where p is the 
mass density. The resulting scatter plots are shown in Figure 13. The straight lines are the 
regression lines. We find very good correlation between the components, valid to above the 
99.9 % confidence level. The slopes are 1.0 (jr), 0.8 (y) and 1.0 (z) for the various compo- 
nents. The positive correlation implies waves propagating against the field. From Figure 12 
we see that average B x > 0 (average = 1.64 nT; sunward) so the waves are propagating 
antisunward. 

A second example is that of 23-31 March 2009. Figures 14 and 15 show similar data 
to that of the previous example. The high-speed stream reached much lower values than in 
the first example. The five-day period 25 - 29 March is analyzed further. Here the average 
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Figure 12 Data at one-minute 
resolution for the period 
29 September -8 October 2008. 
From top to bottom: The proton 
density, temperature, magnetic 
field strength and components, 
the velocity components and the 
AE index. The horizontal bar in 
the last panel gives the eight-day 
interval examined for Alfvenic 
fluctuations. 
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B x is < 0 nT (— 1.5 nT). For waves propagating antisunward, the field and flow compo- 
nents should be anticorrelated, consistent with Figure 15. Again, a reliable anti-correlation 
is found with slopes —0.8 (x), —0.8 (y), —0.9 (z). 

4.2. Shape and Location of Bow Shock and Magnetopause 

We do not propose to examine this topic over the entire 3-year period. To limit our scope, 
we consider Figure 16, which shows in the upper panel daily averages from the OMNI data 
base of the magnetic energy density and in the lower panel daily averages of the kinetic 
energy density. The plot is truncated from below to highlight a feature we are interested 
in, marked by the red arrows. During this about three-month period (27 March -26 June 
2009), corresponding to CR 2081 (part) -2084, we have a clear energy “hole” where both 
energy densities reach their lowest average values. The bottom panel shows the Alfven Mach 
number over the three-year period. 

To characterize the “energy hole” further, Figure 17 shows time profiles of select mag- 
netic field and plasma parameters for the 91 -day period of lowest E mag and £’ klll . From top 
to bottom: the proton number density, bulk speed, total field strength, temperature, dynamic 
pressure and the Alfven Mach number. The data are 5 -min averages. Keeping in mind that 


Springer 


478 


C.J. Farrugia et al. 


Figure 13 The panels show the 
x, y, z components of the 
magnetic field perturbations 
perpendicular to the background 
magnetic field plotted as a 
function of the corresponding 
velocity components multiplied 
by sqrt(pop), where p is the 
mass density. The regression 
lines are shown. R x , y, z. 
represent the respective 
correlation coefficients. 


OMNI 1 min Oct 1 -- Oct 8, 2008 



AA xi (nT) 



AA yi (nT) 



AA 2i (nT) (nT) 


fast co-rotating streams are hot and the slow flows are cold, we have drawn vertical guide- 
lines passing through (most) minima in T in the first five panels (done by eye). These then 
typically occur at stream interfaces (Sis), where the plasma and field are compressed and 
the flow exhibits a sharp positive gradient. This is true even when a slow stream is over- 
taking a still slower stream. Some average values are: N: 5.6 cm -3 ; P (dynamic pressure): 
1.1 nPa (about one-half of its historic average of 2.2 nPa); B : 3.7 nT; M A : 12.0, and B- 
(GSM): — 0.2 ± 2.0 nT (not shown). There is considerable fluctuation about these averages, 
of course, but the values show a solar wind of low kinetic pressure, weak magnetic field and 
small B z component. These features will determine the interaction with the terrestrial mag- 
netosphere during the selected period that we discuss next. Also note the frequent SIRs/CIRs 
and in the later part we have slow-slower SIRs. We now ask: What is the effect of this on 
the shapes and location of the terrestrial magnetopause and bow shock? 

To answer this, we proceeded as follows (see Harris, 2011): We compiled a data set of 
crossings of the magnetopause and bow shock made by Geotail, Cluster 1 and the THEMIS 
B and C spacecraft. These crossings cover both sides of the noon meridian. For Geotail we 
took plasma and magnetic field data from Low-Energy Particle (LEP) experiment (Mukai 
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Figure 14 Similar to Figure 12, 
but for the period 23 - 29 March 
2009. 


OMNI 1 min March 23-31, 2009 
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et al, 1994) and from the Geotail Magnetic Field experiment (Kokubun et al, 1995). For 
Cluster, we used data from the Cluster Ion Spectrometry (CIS) instrument (Reme et al., 
2001) and from the Cluster Magnetic Field Investigation (Balogh et al., 2001). THEMIS 
data are from the Fluxgate Magnetometer (Auster et al. , 2008) and from the Electrostatic 
Analyzer (ESA; McFadden et al., 2008). Interplanetary data were from the Magnetic Field 
Investigation (MFI, Lepping et al., 1995) and the Solar Wind Experiment (SWE, Ogilvie et 
al., 1995) on the Wind spacecraft. This search resulted in a total of 328 unambiguous mag- 
netopause crossings and 271 bow shock crossings. The bulk of the data are in the X-range 
[15, — 20] 7?e- The data were in GSE coordinates. We assumed a paraboloid (hyperboloid) 
of revolution for the magnetopause (bow shock), just as in Fairfield (1971). Each data point 
was then corrected for aberration due to the motion of the Earth around the Sun. To deter- 
mine the solar wind conditions during this period, we used averages over the interplanetary 
parameters shown in Figure 17. When multiple magnetopause or bow shock crossings were 
observed within 30 min, only the average magnetopause position is plotted. The interplan- 
etary conditions correspond to this average time. Note that the three-month period under 
survey contained only one ICME (Richardson and Cane, 2010) and so the effect of ICMEs 
is not present in the result. 


Springer 


480 


C.J. Farrugia et al. 


Figure 15 Similar to Figure 13, 
but for the period 23 - 29 March 


2009. 


OMNI 1 min March 25 — March 29, 2009 
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The IDL fitting routine AMOEBA is used. It computes a multi-dimensional minimization 
of a user-specified function. The expression we minimize is of the same functional form as 
that of Fairfield (1971) and is given by 

0 = y 2 + Axy + Bx 2 + Cy + Dx + E. 

Figure 18 shows in blue the resulting fit and in red the classic Fairfield (1971) result. 
The interplanetary data for Fairfield’s model came from the IMP 1-6 spacecraft in the time 
frame (not continuous) from the end of 1963 to 1968, i.e. in the rising and maximum phase 
of Solar Cycle 20. The fitted coefficients are: 

a = -0.1261, b = 0.0245, c = 0.6717, d = 21.046, e = -248.203. 

Two interesting features may be noted: (1) The stand-off distance of 11.8 Re that we 
obtain is noticeably larger than that of Fairfield’s result, and (2) and the magnetopause flares 
out more in 2007 - 2009. Both are likely due to the low average /Ai yn during this three-month 
period. 

In Figure 19 we compare the fitted magnetopause shape (in light blue) with two model 
magnetopauses, that of Shue et al. (1998; in red) and that of Sibeck, Lopez, Roelof (1991; 
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OMNI 2007-2009 Energy densities 





MONTHS from Jonuory 2007 


Figure 16 Daily averages of the magnetic and kinetic energy densities. Marked by arrows is a three-month 
period where both energy densities are on average the lowest reached in the three years under survey. The 
bottom panel shows the Alfven Mach number over the three-year period. 


in dark blue). The Sibeck, Lopez, Roelof model is based on 1821 crossings of the magne- 
topause. The data are binned by five values of the dynamic pressure. We plotted here the one 
corresponding to the average solar wind dynamic pressure in our chosen period. 

In the Shue et al. (1998) work the influence of both solar P dyn and IMF B z on the magne- 
topause shape are considered. No restriction is imposed on these parameters. Here we used 
the B z -south version of Shue et al. (1998) since the average of this quantity - computed from 
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OMNI 5min March 27 - June 26, 2009 GSM 



Doys after Morch 27, 2009 


Figure 17 Time profiles of select plasma and magnetic field parameters: The proton number density, bulk 
speed, total field strength, temperature, dynamic pressure and Alfven Mach number. Vertical blue guidelines 
are drawn through the minima in the temperature, T . For further details, see text. 


daily values - in our period is slightly negative (—0.2 nT). Shue et al. (1998) postulated a 
new functional form for the magnetopause shape, namely, 

r = r 0 [ 2/(1 +cos(0))]", 

where quantities ro and a are the stand-off distance and the tail flaring angle, respectively. 
Quantity r gives the radial distance measured at an angle 9 with the direction of r. 

The figure shows that both model magnetopauses underestimate the flaring of the mag- 
netopause in our period. Of the two, Sibeck, Lopez, and Roelof (1991) model comes closer 
to reproducing the observed flaring in our three-month period. Both models underestimate 
the stand-off distance by ~ 1 R E . 

The same three-month period in 2009 was employed to determine the bow shock shape 
and location. Figure 20 shows the resulting fit to the data (light blue trace). The fitted coef- 
ficients are: 


a = 0.0252, 6 = 0.4061, c = -0.4576, 

d = 40.374, e = -663.475. 
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Figure 18 The solar minimum 
magnetopause (blue curve) 
plotted in (x, R) coordinates with 
R = (X 2 + F 2 ) 1 / 2 . The red 
curve gives the empirical model 
of Fairfield (1971). 



It is compared with Fairfield’s 1971 model (red). The solar minimum bow shock is less 
flared than Fairfield’s. Interestingly, the subsolar location (bow shock “nose”) is practically 
the same. Two major factors affect the subsolar distance of the bow shock: When Ma in- 
creases, the bow shock moves earthward. At the same time, when the dynamic pressure 
decreases the magnetopause moves sunward and flares out. The stand-off distance of the 
bow shock is proportional to the object size (Farris and Russell, 1994), so that for lower 
Pdyn the bow shock shifts sunward. It may be fortuitous that the two effects cancel each 
other, or that the data are not sufficient to show any remaining small offset. 

In summary, the thinning by 1 R E of the subsolar magnetosheath during this three-month 
period is mainly due to a sunward motion of the magnetopause. 

A small asymmetry is observed near the terminators of the bow shock. The dusk side 
crosses at X = 0 at 26.5 R E while the dawn side crosses the terminator plane at 26.0 R E . 
This asymmetry is small, however, and we shall not consider it further. 

The nose of the minimized function describing the bow shock is at 14.35 R E . With the 
result for the magnetopause, this corresponds to a magnetosheath thickness of 2.6 R E . 
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Figure 19 The solar minimum 
magnetopause (light blue trace) 
compared with the Sibeck, 
Lopez, Roelof (1991; dark blue) 
and the Shue et al. (1998; red) 
model magnetopauses. For 
calculating these model shapes 
the average solar minimum 
values for the dynamic pressure 
and B~ were used. 



X (Re) 


5. Discussion and Conclusions 

We have discussed i) solar wind properties during the extreme solar activity minimum of 
2007-2009, and ii) two major aspects of the magnetospheric response. For the former we 
used data from ST-A and found that the density and magnetic field strength were, compared 
to other solar activity minima, significantly smaller. These values decreased the magnetic 
as well as the kinetic energy densities. The ratio of these energy densities, which is equal 
to the inverse square of the Alfven Mach number, was unusually low. This has an effect on 
the interaction processes with the magnetosphere (see below). For ii), we first considered 
the solar wind electric field and found this to be small (~ 1.3 mV/m). The implication was 
that the response of the magnetosphere, in so far as it depends on interplanetary parame- 
ters, remained linear (no saturation effects). This justified our use of an empirical formula 
- that of Boyle, Reiff, and Hairston (1997) - to find the average cross-polar cap potential 
from interplanetary parameters. The average value turned out to be 37 kV. Finally, we 
also derived empirical shapes of the magnetopause and bow shock using for these cross- 
ings on both sides of noon made by Cluster 1 (on the downside) and THEMIS B and C (on 
the duskside) during a three-month period when both solar wind energy densities reached 


•£) Springer 


Deep Solar Activity Minimum 2007-2009: Solar Wind Properties 


485 


Figure 20 Solar minimum bow 
shock (blue). In red is the 
Fairfield (1971) model bow 
shock. 



their lowest values in the three years under survey. These were augmented by Geotail ob- 
servations. Comparing with Fairfield’s (1971) classic result, there resulted a more flared out 
magnetopause whose stand-off distance (at 11.8 Re) was larger by ~ 1 R E . The empirically 
determined bow shock was less flared out than Fairfield’s. MF1D effects tend to increase the 
flaring of the bow shock (see the Merka et al, 2003, their Figure 10). Thus our result is 
likely due to the fact that the solar wind approached closer to the gas dynamic limit (i.e. Ma 
was higher) than during Fairfield’s period, which corresponded to the rising phase of Solar 
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Cycle 20. The subsolar thickness of the magnetosheath was about 1 R E narrower than that 
of Fairfield. Assuming a normal thickness of 4 R E , this amounts to a shrinkage of ~ 25 %. 

A major protagonist in this study is the solar wind (upstream) Alfven Mach number, 
M a . The importance of this parameter is that the magnetic forces in the MHD momentum 
equation scale as Af~ 2 (Erkaev, 1988; Farrugia et al, 1995). Thus increasing M A leads to a 
magnetosheath in which magnetic forces are weaker. (M A equals infinity is the gas dynamic 
limit.) Typically at 1 AU, quantity M A is ~ 8 — 10. An important exception noted in the 
literature is the low values of order 2-6 achieved in magnetic clouds/ICMEs (Farrugia et al . , 
1995), which are magnetically dominated (low /l) structures. This accounts for the unusual 
magnetosheath properties during their passage at Earth and a concomitant altered interaction 
with the magnetosphere (e.g. Lavraud and Borovsky, 2008). During the long minimum we 
are investigating, M A often reached very high values with an average of 12.9 (see e.g. bottom 
panel of Figure 16). Thus magnetic forces are weak in the solar wind prevailing in these three 
years. This, coupled with the repetitive high speeds reached in what is almost a continued 
stream-stream interaction process, provide favorable conditions for the Kelvin-Helmholtz 
instability to grow, since (a) they weaken the stabilizing magnetic tension forces, and (b) 
enhance a destabilizing factor (the velocity shear) a good fraction of the time, when high- 
speed flows were prominent. Indeed Engebretson et al. (1998), studying the effect of two 
repetitive co-rotating stream interactions over a six-month period, found that auroral Ultra- 
Low Frequency (ULF) wave power in the PC 5 range had a power law dependence on solar 
wind velocity. The ULF power was strongest at the leading edge of high-speed streams. 
They obtained good agreement with a numerical simulation which assumed that the Kelvin- 
Helmholtz (KH) instability was mediating this activity. Arguing by analogy, we hypothesize 
that the KH instability, one of the contributors to viscous-type solar wind-magnetosphere 
interactions, played a major role in transferring solar wind energy and momentum to the 
magnetosphere during these three years. 

We note that in 2007 - 2009 the ICME occurrence rate at 1 AU was low in comparison 
to the previous minimum, and the few ICMEs that were detected had lower field strengths 
and lower speeds (Lepping et al., 201 1). In early 2007, less than one ICME per month was 
observed, though this number increased somewhat at the end of 2008 presumably due to 
deflections of CMEs towards the ecliptic (Kilpua et al., 2009). Due to the low CME/ICME 
speeds, about half of the ICMEs were swept up by CIRs before reaching 1 AU in 2007 and 
early 2008 (Kilpua et al., 2009; further ICME identifications are also provided by Richard- 
son and Cane (2010)). Thus, during 2007 - 2009 the interaction of the solar wind with the 
magnetosphere was clearly not dominated by ICMEs but by CIRs/SIRs. 

In fact, as already mentioned, the interplanetary medium resembled a continued stream- 
stream interaction process. The Alfven waves in the high-speed streams propagating out- 
ward from the Sun are known to be the cause for continued auroral activity, which we also 
observe (Figure 11, bottom panel, and Figures 12, 14). As Tsurutani and Gonzalez (1987) 
pointed out, the B z -fluctuations about zero in the high-speed streams are likely to cause 
reconnection and give rise to this activity (“HILDCAAS”, short for high-intensity, long- 
duration, continuous AE activity). So this reconnection interaction, too, contributed to the 
CPCP. 

In Section 2.2 we also gave results on the energy budget in the interplanetary medium, 
discussing in particular quantity QI, which is the ratio of the magnetic-to-kinetic energy den- 
sities. We can put our results on this in a broader context. The values of QI obtained here can 
be compared with two cases studied in the literature: i) QIs in normal solar wind at 1 AU, 
and ii) QIs in interplanetary coronal mass ejections, which are magnetically dominated struc- 
tures (low proton /3). For i) we consider Osherovich, Fainberg, and Stone (1999) who gave 
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a formula for QI in terms of SSN. This is QI = a + b ■ SSN ( a = 7.6 x 10 -3 ±0.15 x 10 -3 ; 
b = 7.8 x 10 -5 ±0.3 x 10 -5 ). Using SSN numbers for 2007 -2009 from the NOAA/National 
Geophysical Data Center (NGDC) website, we obtain an average value of QI = 0.007, in 
good agreement with our data, particularly for 2007 - 2008 (Figure 8). 

In a series of papers, (Leitner et al., 2009; Leitner, Voros, and Leubner, 2009; Leitner 
and Farrugia, 2010; Leitner, Farrugia, and Voros, 201 1) discussed other aspects of the prob- 
ability distribution functions of QI, which complement the present work and even extend 
it in the detailed analysis of the shape of the distribution functions. They also gave other 
details, such as the distributions of QI in the slow and fast solar winds, separately. These are 
beyond the scope of our work here, but the reader is encouraged to look at these papers for 
more detail. 
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